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The key role of the orbital degree of freedom to understand the magnetic properties of layered 
ruthenates is here discussed. In the G-type antiferromagnetic phase of Ca2Ru04, recent X-ray 
experiments reported the presence of 0.5 hole per site in the d xy orbital, while the d yz and d zx 
orbitals contain 1.5 holes. This unexpected t2 K hole distribution is explained by a novel state with 
orbital ordering (OO), stabilized by a combination of Coulomb interactions and lattice distortions. 
In addition, the rich phase diagram presented here suggests the possibility of large magnetoresistance 
effects, and predicts a new ferromagnetic OO phase in ruthenates. 

PACS numbers: 75.50.Ee, 75.30.Kz, 75.10.-b, 71.10.-w 



The single-layered ruthenate Sr2Ru04 has recently at- 
tracted much attention, both in its experimental and the- 
oretical aspects, since it is an exotic material exhibiting 
triplet superconductivity in the solid state Q. By anal- 
ogy with superfluidity in 3 He, the triplet superconduc- 
tivity is believed to originate from ferromagnetic (FM) 
spin fluctuations enhanced by the Hund coupling [Q. 
Although Knight shift experiments confirm spin triplet 
pairing ||, the symmetry of Cooper pairs remains con- 
troversial, since some experimental results are consistent 
with the existence of a line-node gap Q], in contrast to 
the theoretical candidate with a nodeless p-wave gap. 

In addition to the triplet superconductivity, ruthen- 
ates exhibit complex magnetic properties. When Sr is 
partially substituted by Ca, superconductivity is rapidly 
destroyed and a paramagnetic (PM) metallic phase ap- 
pears, while for Cai.5Sro.5Ru04, a nearly FM metallic 
phase has been suggested. Upon further substitution, the 
system eventually transforms into an AFM insulator || . 
The G-type AFM phase in Ca2Ru04 is characterized as a 
standard Neel state with spin S=l HQ. To understand 
the Neel state observed in experiments, one may con- 
sider the effect of the tetragonal crystal field, leading to 
the splitting between xy and {yz,zx} orbitals, where the 
xy-orbital state is lower in energy than the other levels. 
When the xy-orbital is fully occupied, a simple superex- 
changc interaction at strong Hund coupling can stabilize 
the AFM state. However, recent X-ray absorption spec- 
troscopy studies have shown that 0.5 holes per site exist 
in the xy-orbital, while 1.5 holes are contained in the zx- 
and yz-orbitals ||, suggesting that the above naive pic- 
ture based on crystal field effects is incomplete. This fact 
suggests that the orbital degree of freedom may play a 
more crucial role in the magnetic ordering in ruthenates 
than previously anticipated. 

In this Letter, the multi-orbital Hubbard model tightly 
coupled to lattice distortions is analyzed using numerical 



and mean-field techniques. Since this model includes sev- 
eral degrees of freedom (2-spin and 3-orbital/electron), it 
is quite difficult to study large-size clusters. However, it 
is believed that the essential character of the compet- 
ing states can be captured using a small cluster, through 
the combination of the Lanczos method and relaxational 
techniques. Mean-field approximations complement and 
support the results obtained numerically. An important 
conclusion of our analysis is that the G-type AFM phase 
is stabilized only when both Coulombic and phononic in- 
teractions are taken into account. The existence of a 
novel orbital ordering is crucial to reproduce the peculiar 
hole arrangement observed in experiments S. Another 
interesting consequence of the present study is the possi- 
bility of large magneto-resistance phenomena in ruthen- 
ates, since in our phase diagram the "metallic" FM phase 
is adjacent to the "insulating" AFM state. This two- 
phase competition is at the heart of Colossal Magneto- 
Resistance (CMR) in manganites and, thus, CMR- 
like phenomenon could also exist in ruthenates. 

Let us briefly discuss a localized Ru 4+ ion, in which 
four electrons are contained in the 4c? orbitals. Since the 
crystal field splitting between e g and t2 g orbitals is larger 
than the Hund's rule coupling, the Ru 4+ ion is in the 
low-spin state (S=l): three up and one down (or three 
down and one up) electrons occupy the triply degener- 
ate t2 g manifold, spanned by d xy , d yz , and d zx . These 
t2 g electrons can move to the neighboring Ru sites via 
the oxygen 2p7r orbitals, and, in addition, they are corre- 
lated with each other, and are coupled to the distortion of 
the RuOg octahedron. The above situation is described 
by the Hamiltonian H=Hj s i n +H e i- e i-\-H e i- p i l . The first 
term Akin denotes the hopping of <2g electrons, given by 

-^kin ^ ^■y-y / ^i^ cr ^i+a7 / (T; (1) 

ia77'(T 

where d- lla is the annihilation operator for a t2g-electron 
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with spin a in the 7-orbital at site i (7=xy, yz, and zx), a 
is the vector connecting nearest-neighbor sites, and / is 
the nearest-neighbor hopping amplitude between 7- and 
7'-orbitals along the a-direction via the oxygen 2pir bond, 
given by t* y xy =t* x zx =tj y ^=^^=£=1, and zero oth- 
erwise. Note that only the two-dimensional case is con- 
sidered throughout this paper. 

The second term -ffei-ei denotes the Coulomb interac- 
tions between i2 g electrons, given by 

#ci-ci = U ^2 PhlPhl + u '/ 2 ^2 P'nPh' 

+ j'/2 d l* d L' d h'«' d n'«, ( 2 ) 

where pi ia = d\ 1(J d ill7 , p- n = J^aPnv u is tnc intra- 
orbital Coulomb interaction, U' is the inter-orbital 
Coulomb interaction, J is the inter-orbital exchange in- 
teraction, and J' is the pair-hopping amplitude between 
different orbitals. Note that U=U'+J+J' due to rota- 
tional invariance in orbital space, and the relation J'=J 
holds from the evaluation of Coulomb integrals |J . 

Finally, the third term H \ „y, indicates the electron- 
lattice coupling, expressed as fllPf 

Hq\— ph Q ^ (Q ziPixy ~t~~ QxiPiyz ~t~ QyiPizx) 
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(k/2)J2(Qli + Q 2 m), 



(3) 



where g is the electron-lattice coupling constant, k 
is the spring constant, Q x i=(-l/2)Q3i+(\/3/2)<22i, 
Q yi =(-l/2)Q 3 i-(V3/2)Q 2i , and Q zi =Q 3i . Here Q 2i 
and Q 3 i are the (x 2 —y 2 )- and (3z 2 — r 2 )-mode distortions 
of the RuC-6 octahedron, respectively. The self-trapping 
energy is defined as E p h=g 2 /(2k). Since oxygens are 
shared by adjacent octahedra, the distortions do not oc- 
cur independently, and a cooperative treatment should 
be employed in this problem. A simple way to include 
such an effect is to optimize directly the displacement of 
oxygen ions |J. In practice, considering sites i and i + a, 
the oxygen in between is only allowed to move along the 
a-axis (neglecting buckling and rotations), while apical 
oxygens move freely along the z-direction. 

The above model is believed to provide a good starting 
point to study the electronic properties of ruthenates, but 
it is difficult to solve even approximately. To gain insight 
into this complex system, an unbiased technique should 
be employed first. Thus, in this paper, first a small 2x2 
plaquette cluster is analyzed in detail by using the Lanc- 
zos algorithm for the exact diagonalization, and the re- 
laxation technique to determine the oxygen positions. In 
actual calculations, at each step for the relaxation, the 
electronic portion of the Hamiltonian is exactly diagonal- 
ized for a fixed distortion. Iterations arc repeated until 
the system converges to the global ground state 111) . 
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FIG. 1. (a) Ground state phase diagram for the full model 
with J=3(7'/4. The notation is explained in the text. Other 
choices of J lead to a similar phase diagram (see also Fig. 2). 
(b) Electron configurations in t2 S orbitals for the FM/OO and 
G-AFM/OO states. Solid circles connected by solid curves 
denote the local S=l spin. Thin lines connecting orbitals 
denote the allowed hopping processes, (c) Spin correlation 
for the G-AFM phase, (d) Orbital pattern for FM/OO and 
G-AFM/OO states, (e) Orbital ordering in ruthenates, de- 
duced from the 2x2 result, (f) Mean-field energies for CDW, 
G-AFM, and FM states on the 8x8 lattice and E ph =1.0, with 
lattice distortions assumed from the exact 2x2 results. 

The main result of this paper is summarized in the 
ground state phase diagram shown in Fig. 1(a). There 
are six phases in total, which are divided into two groups. 
One group is composed of phases stemming from the 
U'=0 or i? p h=0 limits. The origin of these phases will 
be addressed later, but first their main characteristics 
are briefly discussed. For i? p h=0, a C-type AFM or- 
bital disordered (OD) phase appears in the region of 
small and intermediate U'. This state is characterized 
by n xy :n yz +n zx =l/2:3/2, where n 7 is the hole number 
per site at the 7-orbital. Hereafter, a shorthand notation 
such as "1/2:3/2" is used to denote the hole configu- 
ration. For large U' , and still E p ^=0, a FM/OD phase 
characterized by 3/4:5/4 is stable, which may correspond 
to Sr2Ru04. On the other hand, for U'=0 and small 
Sph, a "metallic" (M) phase with small lattice distortion 
is observed, while for large -E p h, a charge-density- wave 
(CDW) state characterized by 1:1 was found. In short, 
the G-type AFM phase observed experimentally Q does 
not appear, neither for £ , p h=0 nor for U'=0. 
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Another group includes two phases which are not con- 
nected to either -E p h=0 or U'—O. It is only in this group, 
with both lattice and Coulomb effects being relevant, that 
for intermediate U' the G-type AFM and orbital ordered 
(00) phase with 1/2:3/2 found in experiments [|| is sta- 
bilized. At larger U', a FM/OO phase occurs with the 
same hole arrangement. The electron configurations for 
the FM and AFM phases are summarized in Fig. 1(b). In 
the FM phase, since an S=l spin with S z =+1 is formed 
at each site, the up-spin number is unity at each orbital, 
while the down-spin distribution depends on the orbital. 
In the AFM state, the configuration of double-occupied 
orbitals is the same as in the FM phase, but the single- 
occupied orbital contains 0.5 up- and 0.5-down spins on 
average, since the S=l spin direction fluctuates due to 
the AFM coupling between neighboring S—l spins. How- 
ever, as shown in Fig. 1(c), the spin correlations peak at 
(•7r,7r), indicating the G-AFM structure. Except for the 
spin direction, the charge and orbital configuration in the 
FM/OO phase is the same as in the G-AFM/OO state. 
As shown in Fig. 1(d), a clear ordering pattern includ- 
ing xy, yz, and zx orbitals is suggested for these FM and 
AFM phases Deduced from the present 2x2 re- 

sult, the proposed bulk orbital ordering in ruthenates is 
shown in Fig. 1(e). Note that in this pattern xy-orbitals 
occupy only half the Ru-sites, while zx- (yz-) orbitals 
connect second nearest neighbor xy-orbitals along the x- 
(y-) axis. This structure, quite different from a uniform 
population of the xy-orbitals, is natural from the view- 
point of the kinetic energy gain of i2g electrons, with all 
orbitals contributing to such gain. 

Although the results discussed above are exact, they 
have been obtained using a small-size cluster. Unfor- 
tunately, larger clusters cannot be studied due to se- 
vere limitations of memory and CPU time. However, 
it is possible to consider larger lattices by employing a 
mean- field approximation. In the actual calculations, 
only the diagonal portion of the Coulomb interactions 
in spin and orbital space are included by carrying out 
a standard Hartree-Fock decoupling procedure, working 
with the fixed lattice distortions deduced from the 2x2 
results. In Fig. 1(f), the mean-field energies vs. U' on 
the 8x8 cluster are depicted at -E p h=1.0. It is remark- 
able that with increasing U', the ground state changes 
from CDW, to G-AFM, and finally to FM, in excellent 
agreement with the 2x2 results [|l3]. For intermediate 
values of U', the G-type AFM phase with the hole ar- 
rangement observed in experiments is easily obtained, if 
the lattice distortion to produce the orbital ordering in 
Fig. 1(e) is assumed. For large U', on the other hand, 
the FM phase is stabilized (this phase has not been ob- 
served experimentally yet, and it may correspond to a 
FM insulator). The mean-field generates stable phases 
and their competition are in qualitative agreement with 
the small-cluster results for Ep^O and suggesting 
that the orbital ordering coupled to lattice distortions is 
essential to explain the magnetic structure of ruthenates. 

Now let us discuss intuitively the origin of the complex 



00 pattern observed here. Due to the Hund coupling, a 
local S—l spin is formed by the four electrons, as shown 
in Fig. 1(b). Note that each orbital includes at least one 
spin- majority electron, leading to the cancellation of their 
distortions due to the relation Q x i+Qyi+Qzi=0- Thus, 
the 00 pattern is determined by the choice of orbital oc- 
cupancy of the excess spin-minority electron. Note that 
the energy levels of orbitals occupied and unoccupied by 
this excess electron are — 2E p h and +-E p h, respectively. 
The two-dimensional network composed of xy-orbitals is 
occupied in a bipartite manner by the excess electrons, to 
gain kinetic energy (namely, to avoid the total "freezing" 
of the xy-orbitals). Other excess electrons occupy yz- 
and zx-orbitals, with mobility only along one-dimensional 
(ID) chains, coupled to Q x \ and Q y i distortions, respec- 
tively. These distortions should occur in pairs to cancel 
the in-plane total distortions, since Q x (Q y ) denotes the 
elongation along the x- (y-) direction. Then, half of ID 
chains of yz- and zx-orbitals are occupied, to avoid the 
sites with occupied xy-orbitals. The above considerations 
naturally lead to the orbital pattern Fig. 1(e). 

Another interesting aspect of the phase diagram 
Fig. 1(a) is the appearance of a FM/OD phase adjacent 
to the G-AFM/OO phase, suggesting a competition be- 
tween "metallic" FM and "insulating" AFM states. Since 
this two-phase competition is a key concept to under- 
stand CMR manganites M, the present result suggests 
the possibility of CMR-like effects in ruthenates as well. 
If a small-radius alkaline-earth ion is substituted for Ca or 
Sr, CMR-like phenomena may be observed in Ru-oxides. 
Note that a two-phase competition has already been ob- 
served experimentally in the bilayer ruthenates |lj] . 

To complete the intuitive analysis, the cases of -E p h=0 
and U'=J=0 are discussed separately. First, let us con- 
sider the pure Coulombic model if e i=-ffkin+-ffei-ci us- 
ing the Lanczos method. In Fig. 2(a), the ground state 
phase diagram on the (U'—J, J) plane is depicted by 
comparing the energies for the states classified by the z- 
component of total spin. As discussed in the literature 
of the multi-band Hubbard model jnj, the FM phase 
appears for large J and small U'—J. Although it is dif- 
ficult to analyze transport properties of the FM phase 
in this small cluster calculation, the FM state character- 
ized by 3/4:5/4 is likely metallic, since the orbitals are 
disordered. In fact, there is no clear peak in the orbital 
correlation as shown in Fig. 2(b). For large U'—J, the 
singlet ground state appears, but this is not the Neel, but 
the C-type AFM state, since the spin correlation func- 
tion 5(q) has two peaks at q=(7r, 0) and (0, 7r), not at 
(-7T, 7r), as shown in the inset of Fig. 2(a). Note here that 
this singlet ground state is doubly degenerate, since there 
occurs two types of C-AFM states characterized by the 
peak of 5(q) at q=(0,7r) and (0,7r). This degeneracy 
originates from the geometry of the two kinds of possible 
ID bands composed of yz- and zx-orbitals. Considering 
one of the FM bonds among the S—l spins of the C- 
AFM state, the excess electron in the yz- (zx-) orbital 
can move along the y- (x-) direction, providing a kinetic 
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energy gain, while a superexchange interaction occurs in 
the AFM bonds. Thus, the coexistence of FM and AFM 
bonds in the C-AFM phase is due to the balance of ki- 
netic energy and superexchange interaction. Since the 
orbital correlation function in Fig. 2(c) is composed of 
several orbital patterns of degenerate states, the C-AFM 
phase is also OD. In short, in the 2 x 2-plaquette calcula- 
tion, the singlet ground state of a pure H c \ is C-AFM and 
OD if lattice distortions are not included, in disagreement 
with experiments || . The G- AFM state is not stable in a 
purely Coulombic model due to a Hund coupling energy 
penalization for its mobile electrons, compared with the 
C-AFM state that has better electronic mobility. 
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FIG. 2. (a) Ground state phase diagram for H e \. 
denotes the spin correlation in the C-AFM phase, (b)-(c) 
Fourier transform of orbital correlation function for FM and 
C-AFM phases, respectively, (d) Ground state energy vs. 
Eph for H p h- (e) Average lattice distortion Q vs. E p h- In the 
inset, the electron configuration for CDW state is shown. 

Consider now the pure phononic model _ff p h = -ffkin + 
ffci-ph analyzed by using the relaxation technique. In 
Figs. 2(d) and (e), the ground state energy and the aver- 
age distortion Q are shown as a function of -E p h, where 

Q is defined as Q=(1/4)X) 1 y/Qli + Ofv Note that the 
unit of Q is \ftjk, and k is assumed to be constant here. 
For small E p h, the ground state with small lattice dis- 
tortion is characterized by 1/2:3/2. When E p h is fur- 
ther increased, the distortion abruptly becomes larger, 
and the hole configuration changes to 1:1. In this state, 
to gain lattice energy, two of the three t2g orbitals are 
doubly occupied, leading to the orbital-dependent CDW 
state. The electron configuration is shown in the inset of 
Fig. 2(e), indicating that the orbital population changes 
from site to site. Note that the ground state of ff p h is al- 
ways non-magnetic without the electron correlation that 
forbids double occupancy. Then, the existence of a finite 
Hund coupling, and concomitant Hubbard repulsion, is 



crucial to reproduce experimental results, as long as it is 
supplemented by a lattice distortion. 

Summarizing, to understand the magnetic properties 
of ruthenates, a model with Coulombic and phononic 
interactions has been analyzed by numerical techniques 
and mean-field approximations. It has been found that 
the stabilization of the nontrivial orbital ordering state 
described here is crucial to understand the "1/2:3/2" 
G-type AFM phase observed in experiments. Both 
phononic and Coulombic interactions play an important 
role in this stabilization. In addition, (i) a competition 
between FM/OD and AFM/OO phases suggests the pos- 
sibility of large magnetoresistance effects in ruthenates, 
and (ii) a FM/OO phase is predicted to exist for ruthen- 
ates with large Coulomb interactions. 
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